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Many hybrid models of heavy ion collisions construct the initial state for hydrodynamics from 
transport models. Hydrodynamics requires that the energy-momentum tensor and four-currents 
do not deviate considerably from the equilibrium ideal-fluid form, but the ones constructed from 
transport do not necessarily possess this property. In this work we investigate the space-time picture 
of deviations from equilibrium in Au+Au collisions using a coarse-grained transport approach. 
The collision energy is varied in the range Eiab = 5 — 160A GeV. The sensitivity of deviations 
from equilibrium to collision centrality, and other parameters such as the switching criterion, the 
amount of statistics used to construct the initial state, and the smearing parameter a is investigated. 
For low statistics deviations of from equilibrium are large and dominated by the effect of finite 
sampling. For large statistics the pressure anisotropy plays the most significant role, while the 
off-diagonal components of are small in a large volume during the whole evolution. For all 
considered energies and centralities the pressure anisotropy exhibits a similar feature: there is a 
narrow interval of time, when it rapidly drops in a considerable volume. This allows us to introduce 
an ” isotropization time”, which is found to decrease with energy and slightly increase with centrality. 
The isotropization times are larger than times typically used for initializing hydrodynamics. 


I. INTRODUCTION 

Heavy ion collision experiments at the Relativistic 
Heavy Ion Collider (RHIC) at center of mass energy 
y^SiVAT = 200 GeV per nucleon pair at Brookhaven 
and at the Large Hadron Collider (LHC) at Geneva at 
s/snn = 2.76 TeV have provided evidence of the Quark- 
Gluon Plasma (QGP) existence and proven that it be¬ 
haves almost like a perfect fluid HHS]. Models based 
on ideal hydrodynamics, which requires local thermaliza- 
tion, are very successful in describing the experimental 
data from RHIC and LHC. It was argued that the good 
description of experimental data requires rather small 
initialization times (r < 0.6 fm/c at RHIC) [7]. This 
raised the question about how fast and due to which 
physical mechanisms the QCP can be thermalized. In 
contrast to this requirement of early thermalization, in 
the newer works applying viscous hydrodynamics it is 
claimed that ’’experimental data for bulk quantities can 
be reproduced by hydrodynamic models also for large ini¬ 
tialization times, so no early thermalization assumption 
is needed” [8]. Still, the necessary condition for applying 
hydrodynamics is the vicinity to the thermal equilibrium. 
The recently developed anisotropic hydrodynamics m 
is partially an exception: it can be applied for arbitrar¬ 
ily high pressure anisotropy, but not for arbitrarily high 
off-diagonal components of the energy-momentum ten¬ 
sor. Unlike ideal or viscous hydrodynamics, anisotropic 
hydrodynamics is not yet a well-established ansatz, for a 
discussion see [9]. 
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The question about the degree of non-equilibrium at 
the time of hydrodynamics initialization cannot be an¬ 
swered by hydrodynamics itself. The approach to equi¬ 
librium in heavy ion collisions is actively studied micro¬ 
scopically within two types of models, an overview of 
both being given in m- The first ones apply dualities of 
supersymmetric Yang-Mills gauge theory for calculations 
in the strong coupling limit naE]. Other approaches 
are able to achieve fast thermalization in a weak coupling 
limit, where colliding nuclei are described in the color- 
glass condensate (CCC) framework [141 |T5]. The primary 
effect in CCC leading to fast thermalization is believed 
to be plasma instabilities, such as the chromo-Weibel 
instability [TB]. Both kinds of models predict considerable 
momentum space anisotropies at the times, when hydro¬ 
dynamics is initialized. The aforementioned studies are 
relevant for high collision energies. At intermediate en¬ 
ergies thermalization was studied using transport models 
[30] , where momentum distributions were averaged over a 
(5x5x5) fm^ central cell. In contrast to the work |30] , in 
this paper we study thermalization and isotropization at 
intermediate energies locally in space, paying attention 
not only to pressure anisotropy, but also to off-diagonal 
components of the energy-momentum tensor. Our study 
is particularly relevant for hydrodynamics+transport hy¬ 
brid models that construct initial state from transport 
approaches [IHllIiEIlliaESHlZ]. 

In such models the initial energy-momentum tensor 
and the four-current are obtained from parti¬ 
cles, a procedure we further refer to as fluidization. Af¬ 
terwards they can be unambiguously decomposed into 
= Pleat + and where 

and are the energy-momentum tensor and the four- 
current of conserved charges of an ideal fluid. For a 
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fully consistent fluidization a set of conditions is required. 
First, the hydrodynamics should be applicable, therefore 
the system has to be close to equilibrium. 

ik'^'^ii«i&ii (1) 

m\<^\\jLal\\ ( 2 ) 

For applicability of hydrodynamics it is also required that 
the Knudsen number is small, 

Kn = Ifnicro /L macro « 1 , ( 3 ) 

where Imicro is a microscopic scale related to particle col¬ 
lisions, e.g. mean free path, and Lmacro is a macroscopic 
scale, e.g. the system size. Second, hydrodynamics and 
transport have to be matched: energy-momentum tensor, 
currents and their time derivatives should be the same 
for transport and hydrodynamics at fluidization. Match¬ 
ing time derivatives means that the transport properties 
(viscosities, thermal conductivity, relaxation time, etc), 
chemical rates and the equation of state (EoS) should be 
the same on both sides: 


rpIlV _ rjiflO 

-‘-transport -‘-hydro 
j transport ^hydro ( 3 ) 

same transport properties, chemical rates and EoS ( 6 ) 

Let us call this set of conditions ([T|^ strong consis¬ 
tency conditions. These conditions guarantee that for a 
small time after fluidization results of the hydrodynamic 
calculation are equivalent to the results of the transport 
model. Hydrodynamics in this case is just another lan¬ 
guage for the description of the same phenomena, which 
may have some technical advantages. This is the case for 
the applications of non-relativistic hydrodynamics (see, 
e.g., [28] and references therein): for boundary layers and 
dilute regions, where hydrodynamics is not applicable, a 
transport model is used, otherwise hydrodynamics is ap¬ 
plied, because it is computationally cheaper. In such 
models simulating the whole system with transport and 
comparing to the hybrid model is just a consistency check 
- both approaches should give identical results. Strong 
consistency conditions are a necessary and sufficient con¬ 
dition for that. 

The purpose of applying hydrodynamics in the heavy 
ion collisions simulations is completely different from the 
non-relativistic applications, where the hydrodynamics 
is just a matter of the simulation speed-up. In heavy 
ion collisions switching to hydrodynamics allows to in¬ 
vestigate the effects of thermalization, equation of state 
and transport properties of the hot and dense hadronic 
or QGP matter. This implies instantaneous thermaliza¬ 
tion, rapid change of EoS and transport properties at 
fluidization. Therefore, conditions (§ and § are not 
required, but conservation laws represented by equations 
ii should be fulfilled and conditions ([T][^ remain as 
a practical requirement to solve the hydrodynamic equa¬ 
tions numerically. We further call conditions ( [l|2|4| ^ 
the weak consistency conditions. 


The fluidization is typically performed either at a con¬ 
stant proper time hypersurface r = const or at a constant 
center of mass frame time hyper surface tcM = const. 
The constant is often chosen according to the geometri¬ 
cal criterion - the time, when nuclei geometrically pass 
through each other: tcM = = 2R{Eiab/‘imN)~^^‘^, 

where R is radius of the nucleus, /3 is velocity, 7 = 
(1-/3^)“^/^, Eiab is laboratory frame kinetic energy per 
nucleon, and ruAr is nucleon mass. This time is taken 
to be the same for all collision centralities. It was never 
systematically verified, if consistency (strong or weak) is 
fulfilled at fluidization. 

In this contribution we want to close this gap and an¬ 
swer the following questions: 


• Which conditions of weak consistency are fulfilled 
at the time of geometrical overlap? 

• What is a better fluidization criterion? Should it 
depend on centrality? 

• How can one perform fluidization consistently in 
the weak sense? 


To answer these questions we analyse the outcome of 
a hadron cascade and determine, where and when condi¬ 
tions are fulfilled. 


The article is organized as follows. To set the stage, the 
way fluidization is performed in different models is dis¬ 
cussed in Section [IT] In Section [TTI| we derive the expres¬ 
sions to quantify deviations of and from the ideal 
fluid. The methodology of our work is described in Sec¬ 
tion [IV| Results for quantifying the deviations are split 
into two sections: in Section jV] we consider effects of nui¬ 
sance parameters: statistics, smearing and grid spacing, 
while in Section [Vl| we highlight the dependence on colli¬ 
sion energy and centrality. Our criticism to the smearing 
procedure in the fluidization of the existing models and a 
suggestion for its improvement are given in the Appendix. 


II. OVERVIEW OF FLUIDIZATION IN 
CURRENT HYBRID MODELS 

Let us start by summarizing the current approaches 
that are employed in hybrid models (see Table [^. All the 
shown approaches need to obtain the ideal fluid part of 
and from discrete degrees of freedom (hadrons, 
partons, strings). The viscous corrections are neglected 
in all models, even if viscous hydrodynamics is applied 
for the evolution. The only exception is a recent work by 
Liu et al. HU, where the initial stage for viscous hydro¬ 
dynamics is constructed from free streaming partons and 
viscous corrections are explicitly included. 

The energy-momentum tensor and four-currents 
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are constructed as 

TZti'r)= ( 7 ) 

jLtir)= 

Here K is a smearing kernel, which is often taken as 
a Gaussian in various coordinates. In the Appendix we 
argue, that K{r)d^r should be Lorentz scalar, show that 
none of the models fully meets this requirement, and de¬ 
rive a simple smearing kernel, which can be employed 
instead. Note that changing to this Lorentz-invariant 
kernel leaves the final results rather unaffected. 

There are two ways in the literature to match the ob¬ 
tained and to ideal hydrodynamics. The first one 
is to use only assuming that they have ideal 

fluid form = {e + p)u^^u'' - pg^’', = nu^^, 

and adding the equation of state (EoS) p = p{e^n). The 
following system of equations is then solved (usually it¬ 
eratively, for details see EZl): 

rT00 = (e+p)72-p 
T^ = (e + p) 72 u (8) 

3 b = 

\p=PEos{n,e) 

The advantage of this method is that it conserves en¬ 
ergy and momentum. However, this method supports 
switching only to ideal fluid keeping viscous cor¬ 

rections is hardly possible. Even though the switching 
method conserves energy and momentum, one of the 
models, which employs it, violates conservation laws EH, 
because in m the whole is multiplied by a free pa¬ 
rameter which is then fixed by experimental multi¬ 
plicities. 

Another way is to determine the energy density e and 
the collective velocity by solving the eigenvalue prob¬ 
lem 

= eu^ , (9) 

using the fact that is a timelike eigenvector of 
and satisfies = 1. Then the density n is computed 

as n = j^u^. Only after that the pressure is determined 
from EoS. Note that this way is not equivalent to the pre¬ 
vious one: here the collective velocity does not depend 
on the equation of state. This method conserves energy 
and momentum only if the viscous corrections are kept. 
If they are neglected (as in [2U|22l[2S]) then conserva¬ 
tion laws are violated. Eor a simple example assume that 
= 7(l,0,0,i;). In this case, the energy density in the 
computational frame is Ccomp = 7^(^ + where 

rf can be split into the ideal fluid pressure and a vis¬ 
cous correction. If the correction is neglected, energy 
conservation is violated. 

The switching criteria are rather similar in most of the 
approaches: either constant proper time r or constant 


time in the center of mass frame tcM^ the value being 
determined by the geometrical criterion. It provides the 
earliest time, when equilibration could be in principle 
possible. In [20] the fluidization time tfi has a better 
physical motivation: tfi is chosen such that the entropy 
per baryon does not change any more at t > tfi. 

III. EXPRESSIONS FOR VERIFICATION OF 
WEAK CONSISTENCY 

Let us rewrite the conditions of weak consistency with 
hydrodynamics for and in a way convenient for 
numerical computation. General expressions for and 
in viscous hydrodynamics (Landau picture) are the 
following: 

- A^^(Po + H) + (10) 

, 

where H is the bulk pressure, is the shear stress ten¬ 
sor, no is the conserved quantum number density and 
is the diffusion current. Viscous corrections to ideal 
hydrodynamical and are supposed to be small: 

H <C Po (12) 

Ike'll « no (13) 

From Eqns. [^one obtains 

^ _ gp) (14) 

o 

Po + n = -iA^,r^" (15) 

g'' = AOr (16) 

One can see that in the Landau rest frame = 
diag(l, 0, 0, 0), = 0, and = 0. The non-zero com¬ 

ponents are written as follows: 

Po + n = i(ry+ Tf+ Tf) (17) 

7 = n^'-(Po + n)5*^' (18) 

= “Jl (19) 

Let us note that tensor and vector norms are frame- 
independent, so the consistency conditions for viscous 
hydrodynamics can be formulated in any frame. In Eqn. 

one can substitute ||T^^|| by its largest component 
in the Landau frame: eo- Then Eqn. 0 will turn into 

I |Tr - diag(eo, P',P', P') 11 « eo , (20) 

where P' denotes \{Tj} + + T|^) = Pq + H. The 

physical meaning of this equation is that the diagonal 
components of in the Landau rest frame do not de¬ 
viate much from P' and simultaneously off-diagonal com¬ 
ponents are small compared to eo- The condition for q^ 
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Model 

Initial condition 

Hydro 

Switching 

criterion 

Smearing 

kernel 

Getting 

rjijJiU 

^ideal 

UrQMD 
hybrid [H] 

UrQMD 

cascade 

ideal 3+ID, 
SHASTA 

tcM [fm/c] = 
max{2R^ 

Gaussian 

z-contracted 

o 

o 

Skokov-Toneev 
hybrid [El 

Quark-Gluon- 
String-Model 

ideal 3+ID, 
SHASTA 

tcM such 

that S/Qb — const 

not 

mentioned 

o 

o 

EPOS [n] 

Strings (Regge- 
Gribov model) 

ideal 3+ID 

T 

Gaussian 

z-contracted 

Landau frame 

NeXSPheRIO 
hybrid [22l [23] 

Strings (Regge- 
Gribov model) 

ideal 3+lD, 

SPH 

T = 1 fm [^ 

Gaussian in 

X, y, rr] 

Landau frame 

Gale et al [25] 

IP-glasma 

viscous 3+ID, 
MUSIG 

T = 0.2 fm/c 
(v/W = 2.76 TeV) 

not 

mentioned 

Landau frame 

Karpenko 
hybrid [26] 

UrQMD 

cascade 

viscous 3+ID 

'T'geom 

Gaussian with 
a± and arj 

rpnQ ^-0 

Pang et al 
hybrid [27] 

AMPT 

ideal 3+ID, 
SHASTA 

r 

Gaussian with 
a± and arj 

O 

o 


Table L Fluidization features in different hybrid approaches. Each of these models, including those using viscous hydrodynamics, 
neglects viscous corrections at fluidization. 


is rewritten as 


|Ty -Tf+ry -Tf 1/3 < |Ty -Tf 1/3+|ry -rf 1/3 

and Co is substituted by by P'. In this way a set of criteria 

( 21 ) 

is obtained that we use for numerical calculations. 


Here the physical meaning is that relative velocity be¬ 
tween Landau and Eckart frames should be small. To 
rewrite H ^ Pq one has to add an equation of state 
^0 = lo the system. Then one obtains 

P'/PEosi^o^ Jl) — 1 1. Consequently, whether the 

tensor is suitable for fluid dynamics or not is also 
defined by the equation of state from the fluid dynamics 
itself. The same can be consistent with viscous hy¬ 
drodynamics with some equation of state, and may fail 
when the equation of state is changed. Therefore, we will 
not study the smallness of bulk corrections further, but 
leave this for a future study. 

The conditions for smallness of the shear stress ten¬ 
sor can be split into two: pressure isotropy and small¬ 
ness of off-diagonal elements. One has to note that the 
Landau frame is defined only up to an arbitrary rota¬ 
tion. Locally one can always choose such coordinates 
that = T‘1^ = = 0. However, our coordinates 

are the global coordinates of the computational frame 
and therefore non-diagonal components of the are in 
general non-zero. Therefore, 


X = 


IT/I - Tf I + |Tf - Tf I + |Tf - Tj 




rpll 

-^L 


Y = 


Tf 


Tf 


3(|T22| + |Tf I + \Tl^\) 


Til 


■TF 




<1 (24) 

< 1 


VlE 


^ ry + Tf+ rf _ 

“ SpEosieoJl) 


In the following X is referred to as pressure anisotropy 
and Y as off-diagonality. Please note that due to the 
inequality \a — b\ < |a| + |6| it is always fulfilled that 
X < 2. For ideal fluid dynamics X = 0. For Y let us 
remark that 


rj-i\2 


pO 


^ 1 Pi +P2 
- 2 ^ pO 


t(T“+T^y (25) 


Interchanging indices and substituting this into the defi¬ 
nition of Y one gets Y < 3. For an ideal fluid Y = 0. 


|Ty - P'\ + \Tl^ - P'\ + |Tf - P'l < eo (22) 

|Ty| + |Tf| + |Ty|«eo (23) 

To make these conditions stricter, every term is substi¬ 
tuted by the right hand side of the inequality — P'| = 


IV. METHODOLOGY 

Our calculation is based on the hadronic transport ap¬ 
proach - Ultrarelativistic Quantum Molecular Dynamics 
(UrQMD 3.4) [29]. The degrees of freedom in UrQMD 













5 


are hadrons, resonances up to a mass of 2.2 GeV and 
strings and the implemented processes include binary 
elastic and inelastic scatterings which mainly proceed via 
resonance formation and decays or string excitation and 
fragmentation at higher collision energies. The UrQMD 
particles move along classical trajectories and scatter ac¬ 
cording to their free-particle cross-sections. In our stud¬ 
ies there are no long range potentials and particle trajec¬ 
tories between collisions are always straight lines. Using 
UrQMD we simulate An + An collisions at laboratory 
frame energies ^lab = 5, 10, 20, 40, 80 and 160A GeV. 

The general procedure for our calculations is: 


1. Generate many UrQMD events and coarse-grain 
them using a 2-MD space-time grid. The space 
dimensions are chosen to be the event plane xz. 
We choose a 2+lD grid and not 3+lD, because 
observing the behaviour of some quantity on 2D 
surface versus time is much easier and informative 
for a human than observing quantity on a 3D grid. 
Additionally, in central collisions due to symme¬ 
try the event plane completely characterises all the 
volume. The center of mass frame is used as the 
computational frame in the simulation. In all the 
following text the time is measured in the center 
of mass frame. By the UrQMD convention t = 0 
is the moment when the contracted spheres of the 
nuclear radius first touch each other in a central 
collision. 

2. Use particles from the generated events to con¬ 
struct the energy-momentum tensor x, z) lo¬ 

cally for each grid cell. To compute we take 
only participants into account, i.e. the particles 
that took part in at least one collision. To con¬ 
struct Gaussian smearing is employed, for de¬ 
tails see Appendix. We argue that the construc¬ 
tion of in the existing models is not a Lorentz- 
invariant procedure due to the smearing kernels. A 
simple kernel is suggested, which satisfies the nec¬ 
essary Lorentz-transformation requirements. How¬ 
ever, the introduction of the new kernel appeared 
to be only a matter of physical rigour: switching 
standard kernel to our kernel did not produce no¬ 
ticeable changes in the results. 

3. Transform the constructed energy-momentum ten¬ 
sor in each grid cell to the Landau rest frame. 

4. Verify the weak consistency conditions (see Eqn. 


24) locally in time and space. We check conditions 


for smallness of pressure anisotropy, off-diagonality, 
Eckart velocity relative to Landau, but the con¬ 
dition for smallness of bulk pressure compared to 
pressure is left for future studies. 


In this procedure we follow exactly the fluidization in 
hybrid models, the only addition is that the weak consis¬ 
tency checks are performed. 



number of pions N 


number of pions N 


Figure 1. Pressure anisotropy X (left) and off-diagonality Y 
(right) of for particles sampled according to thermal dis¬ 
tributions. The effect of statistics on the deviation of energy- 
momentum tensor from the ideal fluid one is demonstrated. 


V. SENSITIVITY TO STATISTICS, GRID 
SPACING AND SMEARING 

The deviations of in the simulation from the ideal 
fluid can have two distinct reasons. The first one 

is the deviation of the distribution function /(r,p) in 
the transport approach from equilibrium, this reason is 
referred to as physical. The second reason is statistical: 
due to finite number of particles in the simulation, the 
distribution function is not sampled exactly. 

Let us consider two energy-momentum tensors: calcu¬ 
lated from particles and a ’’true” 

^part(r) = J^ E ( 26 ) 

events i 

= / ^f(r,p)d^P (27) 

Here Ney is the number of events and K{r — Vi^pi) is a 

smearing kernel. In the limit of N^v oo 

X (28) 

events i 

j d^pdi^r f{r',p)K{r — r') 

Eor the case of a Gaussian kernel for cr ^ 0, the ker¬ 
nel K S{r — r') and one retains the ’’true” To 

combine these two limits (a ^ 0, Ney oc) one has to 
keep enough particles within a volume of size . Conse¬ 
quently, to obtain the ’’true” in the simulation, one 
has to take the limit {p is density) 

(j —y 0, Neypcr^ —y oc (^9) 

This creates practical limitations for finding the ’’true” 

in simulations: decreasing a 10 times demands in¬ 
creasing the statistics 1000 times! One can also see that 
regions with lower density are more demanding with re¬ 
spect to statistics. To get some insights into the effect of 
statistics, we performed an auxiliary simulation: N pions 
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Figure 2. Event plane area percentage, where pressure 
isotropy X (left) or off-diagonality Y (right) does not exceed 
0.3 versus time for different number of events Nev used to 
construct . Number of events Nev — 1 corresponds to 
event-by event case. 




Figure 3. Event plane area percentage, where pressure 
isotropy X does not exceed 0.3. Au+Au versus time. E — 
80A GeV, centrality 6 = 6 fm, number of events Nev — 1000. 
Gaussian smearing a (right) and grid spacing Ax = Xz (left) 
are varied to study sensitivity of results to them. 


are generated, their momenta being sampled from a ther¬ 
mal distribution with an ad-hoc temperature of T = 0.2 
GeV, then is computed, pressure isotropy X and 

off-diagonality Y of the energy-momentum tensor from 
Eq. [^are calculated. We varied number of pions N and 
plotted X and Y versus N. The results can be seen in 
Fig-B For every point the simulation was repeated 100 
times and the standard deviation is displayed as an error. 

Fig. [T] can be used to specify the number of events 
needed to reach a good enough approximation to the 
’’true” For example, for Y = 0.3 as an acceptable 

level, the condition of Eqn. 29 becomes Neypa^ > 100. 
From Fig. [^one can also see that the off-diagonality Y 
is more sensitive to statistics than the pressure isotropy 
X. 


In the previous paragraph we considered the effect of 
statistics itself rather as an obstacle to get the physical 
’’true” However, recently event-by-event simula¬ 

tions gained popularity, where the initial state for hydro¬ 
dynamics is intentionally constructed from a small num¬ 
ber of events to include the fluctuations. Let us see, how 
the number of events influences deviations of from 
the ideal form in heavy ion collisions. We compute 
locally on every point of the grid, and as a general charac¬ 
teristic we choose the percentage of the event-plane area, 
where X < 0.3 {Y < 0.3). To define the total area numer¬ 
ically, only grid cells, where pressure p > 10“^ GeV/fm^ 
are taken into account. For this example Au-fAu colli¬ 
sions at E = 80A GeV with the impact parameter 6 = 6 
fm are considered. The smearing a is 0.8 fm. Results are 
depicted in Fig. 

One can see that for this given a 1000 events are 
enough for X to saturate, so the line for N^v = 1000 
represents results for the physical pressure isotropy, i.e. 
due to deviation of f{r,p) in the transport from equilib¬ 
rium. For Y at Ney = 10000 almost all the event area 
has small off-diagonality, which means that the physical 
off-diagonality is small. For event-by-event simulations 
deviations of from ideal fluid are dominated by sta¬ 
tistical effects. 



Figure 4. Isotropization time (such time that more than 50% 
of event plane area have pressure isotropy X < 0.3) versus a. 


In addition to the effect of statistics we investigate the 
effect of other nuisance parameters, i.e. the grid spacing 
and the Gaussian smearing, on our results. According 
to Fig. the grid spacing does not influence the results 
if taken sufficiently small. This is expected, because the 
grid does not participate in the simulation or in the calcu¬ 
lation of it only determines the resolution of the 
output. The only effect of grid spacing is on the precision 
of the area calculation by counting cells, here resolution 
of the output matters. At early times, it makes some 
difference, because the total area is small. That is why 
further we take Ax = Az = 0.6 fm. At the same time 
the Gaussian a influences the results very significantly, 
as it can be seen from Fig. The effect of Gaussian 
a is twofold: on the one hand larger a means effective 
increase of statistics. On the other hand, if the pressure 
anisotropy is large at some space point due to physics, 
the Gaussian smearing will spread this asymmetry in a 
1-2 a radius. 

To characterize the influence of a in a simpler way, we 
plot tiso versus a, where Uso is the earliest time when at 
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1^11_ rp 22 |_|_|^22_ rp 33 \\\^33 _^111 

Figure 5. Space-time evolution of pressure anisotropy X = ^ Liif^^ 22 fJ 33 ^ ^ (see color scale above the Fig.) for 

collision energy in lab frame F' = 80^ GeV, centrality 6 = 6 fm. If the value of X exceeds color map maximum, it is marked 
with the same color as maximum. Solid lines mark the positions of the nuclei, if they wouldn’t interact. 


least 50% of the area have X < 0.3. We refer to this time 
as isotropization time as it is further described in the fol¬ 
lowing Section. In Fig.j^this dependence is displayed, the 
isotropization time is monotonously growing with a and 
is approaching to the geometrical time for a ^ 0. Taking 
the limit a ^ 0 is computationally challenging, because 
one has to increase statistics as as we have shown 
previously. Instead we choose a reasonable a = 0.8 fm 
and assign systematic errors to our results, correspond¬ 
ing to changing a in the range (0.6 — 1.0) fm. Another 
justification for such a treatment is that none of the ex¬ 
isting models attempts to consider the ’’physical” limit 
of cr ^ 0, Nevpcf^ oc, all the models use some fixed 
sigma instead. 

VI. RESULTS AND DISCUSSION 

While in the previous section we studied the effects 
of nuisance parameters on the energy-momentum ten¬ 
sor generated from particles, here we consider the de¬ 
pendence on physical parameters: collision energy and 
centrality. All the following figures are shown for grid 
spacing Ax = Az = 0.6 fm, Gaussian smearing cr = 0.8 
fm, and number of events = 1000. 

The pressure anisotropy X satisfies the following 
properties: X < 2 for any tensor and X = 0 if and only 
if T}} = as it is for ideal fluid For viscous 

hydrodynamics it is necessary that X ^ 1. We consider 
Xcrit = 0-3 as a limiting value, when viscous hydrody¬ 
namics is still applicable. Changing Xcrit to 0.4 leaves 
us with qualitatively the same results and conclusions. 
Fig. m gives a qualitative impression of the space-time 
evolution of the pressure anisotropy. Even though the 
figure shows a particular energy of Eiab = 80A GeV and 
centrality b = 6 fm, some features are universal for all 
energies and centralities that we considered: 

• On the borders of the expanding system the 
anisotropy is always high, thus these regions are 
never consistent with viscous hydrodynamics. 



Figure 6. Example of pressure anisotropy behaviour versus 
time. Gollision energy in lab frame E = 80A GeV, centrality 
6 = 6 fm. 

• After some moment of time a relatively isotropized 
central region rapidly expands and never disap¬ 
pears completely during the time evolution 

To make quantitative statements let us consider the 
evolution of the pressure anisotropy at several points 
along the x axis (z = 0) versus time. This is shown in 
Fig.m One can see that in the beginning the anisotropy 
is almost maximal, then it rapidly decreases and never 
rises too much again. As we have shown in the previous 
section, Nev = 1000 is enough to suppress the anisotropy 
due to statistics, so the behaviour of X in Eig.j^is domi¬ 
nated by physics. However, the impact of statistical fluc¬ 
tuations can be observed already at t > 5 fm/c: X starts 
to fluctuate in space and time. One can see this both in 
Eig. [^at later times and in Eig. 

After looking at the space-time evolution of the 
anisotropy in 2D as in Eig. for different energies, one 
gets the impression that there is a special moment Uso 
for each energy and centrality, before which pressures 
are highly anisotropic in the whole event plane and after 
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Figure 7. Percentage of area in the event plane, where 
pressure anisotropy X < 0.3, for Au+Au collision energies 
Eiab = 5 - 160A GeV. 


which there emerges a considerable isotropic region. To 
quantify this feeling we consider the ratio of area, where 
X < 0.3 to the total area versus time. From Fig. [^one 
can see that there is indeed a steep rize of isotropized 
area at some moment for every considered energy and 
centrality. Let us define the isotropization time tiso such 
that more than 50% of the area has X < 0.3 at t = tiso- 
The behaviour of this isotropization time versus energy 
and centrality is compared to the geometrical criterion 
in Fig. One can see that the isotropization time de¬ 
creases with energy, but remains larger than the geo¬ 
metrical criterion for all energies except 5A GeV. It is 
interesting to note that the isotropization time differs 
with centrality: for larger impact parameters b it slightly 
increases. One might assume that this has a pure geo¬ 
metrical reason: to = 0 is chosen in UrQMD as a moment 
when nuclei touch each other in a central collision. How¬ 
ever, for peripheral collisions the nuclei will only touch 
ai to{b) = ^{1 — a/1 — {b/2RY). In Fig. one can see 
that this naive expectation yields the rignt trend: tiso 
rises with centrality and the rise is smaller for higher en¬ 
ergies. However, quantitatively it overestimates tiso for 
large impact parameters. 

Let us compare our findings to the study by Bravina 
et al. m, where one central cell of (5 x 5 x 5) fm^ was 
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Figure 8. Isotropization time tiso (see definition in the text) 
versus energy and centrality. 



Figure 9. Isotropization time tiso (see definition in the 
text) versus centrality. Dotted lines are naive expectation 
from the collision geometry: tisoip) — tisoib = 0 ) + X(i _ 

- {h/2RY). 


chosen to study the pressure anisotropy of the energy- 
momentum tensor in Au+Au collisions. An isotropiza¬ 
tion time was defined, and it did not change significantly 
after zooming the central cell to (1 x 1 x 1) fm^. This 
allowed for the conclusion that isotropization happens 
rapidly in a large volume. The isotropization time deter¬ 
mined in this study also decreases with collision energy. 
All of these results are confirmed in the present study. 
However, our isotropization times are smaller than the 
ones obtained by Bravina et al. There are two possible 
reasons for that. First, we do not count spectators for 
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in lab frame E — A GeV, centrality 5 = 6 fm. If the value of Y exceeds color map maximum, it is marked with the same 
color as maximum. Solid lines mark the positions of the nuclei, if they wouldn’t interact. 



Figure 11. Space-time evolution of relative velocity between Landau and Eckart frames vle — + Oi)^ + (see 

color scale above the Fig.) for collision energy in lab frame E — 80A GeV, centrality 5 = 6 fm. If the value oivLE exceeds color 
map maximum, it is marked with the same color as maximum. Solid lines mark the positions of the nuclei, if they wouldn’t 
interact. 


constructing pressure anisotropies. Second, the criterion 
for isotropization we use is less strict: while we only re¬ 
quire at least 50% of event plane area to have X < 0.3, 
the central cell study demands VzjVx — 1 < 0.1 in the 
whole cell, which corresponds to X < 0.065. 

The off-diagonality Y = char¬ 

acterizes the size of the off-diagonal components of the 
stress tensor compared to the pressure. In ideal hydro¬ 
dynamics Y = 0. For the applicability of viscous hy¬ 
drodynamics it is necessary that Y We consider 

^crit — 0.3 as a value, after which viscous hydrodynam¬ 
ics is hardly applicable. An example for the space-time 
evolution of Y is given in Fig. From this Figure it 
can be seen that in the central region Y is always small. 
On the boundaries Y is typically large due to statistical 
effects. A quantitative study similar to the study of the 
pressure isotropy X shows that for all considered ener¬ 
gies and centralities more than 80% of the event plane 
area have Y <0.3 for the whole time of the evolution. 

The relative velocity between Landau and 
Eckart frames for baryon charge vee is shown in 
Fig. At high enough statistics the relative velocity 
between Eckart and Landau frames is not an important 


factor. It is significant only on the borders of the system, 
where the density is small and statistical effects play a 
role. But in all the rest of the volume, for all the consid¬ 
ered time evolution it remains small. 

The effect of momentum-space cuts 

Previously, we included all participants into cal¬ 
culation. However, it is generally believed that soft par¬ 
ticles at midrapidity thermalize faster, therefore it might 
be insightful to impose cuts in momentum space, if one 
wants to obtain a more isotropic On the other hand, 
applying tranverse momentum pt and rapidity y cuts on 
a perfectly symmetric distribution results in an asym¬ 
metry. In addition, these cuts decrease the statistics, 
which leads to an increase of the anisotropy. To study 
the effect of the kinematic cuts on the space distribution 
of the pressure anisotropy X, we constructed only 
from particles with pt < 2 GeV and \y\ < 1. The ef¬ 
fect of cuts on X over space is shown in Fig. We 
made sure that the statistical effect does not play a sig¬ 
nificant role by checking that results do not change after 
increasing the number of events from Nev = 10^ to 10^. 
The large anisotropy at very early times decreases after 
imposing cuts. But the anisotropy at later times is strik- 
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Figure 12. The effect of pT < 2 GeV and \y\ < 1 
cuts on the space distribution of pressure anisotropy X — 

(see color scale above the Fig.) 


^ 22 , 


^33 


for collision energy in lab frame = 80^ GeV, centrality 
6 = 6 fm. If the value of X exceeds color map maximum, it 
is marked with the same color as maximum. Solid lines mark 
the positions of the nuclei, if they wouldn’t interact. 


ingly larger with cuts, first of all in the regions behind 
the nuclei. While without cuts we were able to introduce 
an ”isotropization time”, when X is smaller than 0.3 at 
more than 50% of the area, with the cuts the anisotropy 
is so high all over the space that the isotropization time 
cannot be introduced anymore. 


VII. CONCLUSIONS 


equivalently) is bound to deviate strongly from the ideal 
fluid form. The off-diagonality appears to be mostly pro¬ 
duced by this statistics effect. For large statistics Y tends 
to be small in all the collision region at all times. The 
pressure anisotropy does not vanish at large statistics, it 
is a physical effect related to the anisotropy of the un¬ 
derlying distribution function f{r,p). As a consequence, 
the initial state from UrQMD with enough statistics is 
suitable for anisotropic hydrodynamics. 

Unfortunately, all the results depend on the smear¬ 
ing parameter a. With larger a isotropization is reached 
later, but the degree of isotropization is higher. From a 
practical point of view that means that selecting large 
a one has to take larger fluidization time. This is in 
agreement with conclusion of m that larger fluidization 
time should be taken for larger a to obtain the same 
pion yield. However, strictly speaking, the physical limit 
is cr ^ 0, a^pNey oo. We found that at small a the 
isotropization time approaches to the time of geometrical 
criterion tgeom = 2i?(F^za6/2mAr)~^^^, so the geometrical 
criterion is partly justified, but only in the above men¬ 
tioned limit. In the existing models the smearing param¬ 
eters and statistics are such that at tgeom Ihe anisotropies 
are very high. 

For the pressure anisotropy X it was observed, that 
it exhibits a similar pattern for all the considered col¬ 
lision energies and centralities: there is a narrow inter¬ 
val of time, when it rapidly drops in a considerable vol¬ 
ume. This feature allowed us to introduce and study the 
isotropization time tiso- It decreases with collision en¬ 
ergy, following the same trend as the geometrical overlap 
time. Based on this finding a new fluidization criterion 
is suggested: Uso = tgeom{E) -b Ato(cr), where Ato de¬ 
pends only on the Gaussian smearing a and can be de¬ 
termined from Fig. We observed a slight dependence 
of isotropization time on centrality: it increases with im¬ 
pact parameter, but the slope of this increase becomes 
smaller and smaller for higher energies. This behaviour 
has a simple geometrical interpretation. 

Testing the smallness of the bulk pressure compared 
to pressure left beyond the scope of this article. As we 
have mentioned, the extracted bulk pressure depends on 
the equation of state, therefore this check has to be per¬ 
formed with a set of EoS used in the literature. Since 
in this article we restricted ourselves to testing the weak 
consistency, the next possible step is to test locally the 
fulfilment of the strong consistency. 


We have studied the energy-momentum tensor from 
the UrQMD transport approach locally in space and time 
and quantified its deviation from the ideal fluid form with 
two numbers: the pressure anisotropy X and the off- 
diagonality Y. First of all, we have shown that X and Y 
depend on the number of UrQMD events Nev used to con¬ 
struct Low statistics implies large deviations, even 
if the underlying distribution function is completely ther¬ 
mal and isotropic. An initial state constructed from less 
than few hundred events (or few hundred testparticles 
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VIII. APPENDIX: AND FROM 

PARTICLES 

The general formula for calculating and from 
particles is 

Ttniti.’’’) = J (30) 

jLtir) = J ^fir,p)d^p 

Here / is the 1-particle distribution function. For a dis¬ 
crete set of particles it is 

fir,p) = ^ S^{p - Ppart)S^{r - rpart) ( 31 ) 

part 

For numerical calculations the spatial dependence of the 
distribution function needs to be smeared in some way, 
therefore the delta-function is substituted by some smear¬ 
ing kernel K(r — Vpart)- The smearing kernel should 
satisfy three conditions: K{r)d^r should be Lorentz 
scalar, it should be normalized as j K{r)d^r = and it 
should approach the delta-function for smearing param¬ 
eters approaching to zero. Surprisingly, all current pop¬ 
ular choices of K are such that K{r)d^r is not a Lorentz 
scalar: 

nil • ^ /l/AU,rGAU 

• Cell averaging: K[r) = < . Here 

I 0, otherwise 

K{r)d^r is not a Lorentz scalar, since volume V is 
not contracted. 

• Caussian averaging with Lorentz contraction in z 

direction: K(r) = Nexp(^— ^ ^ is behav¬ 

ing properly only under boosts in z direction. 

• Caussian in x, y, rj coordinates also behaves prop¬ 
erly only under boosts in z direction. 


Further we derive a simple, not too computationally de¬ 
manding kernel, which satisfies all aforementioned con¬ 
ditions, in particular K{r)d^r being a Lorentz scalar. 
In the rest frame of the particle we take a Gaussian: 
K{rrest) = If Ax'^ are the 

coordinates in the computational frame and Ax^ are the 
coordinates in the rest frame then Ax^ = K^Ax'^ ^ where 
A is the matrix of the Lorentz transformation. We want 
to consider a Gaussian at fixed time in the computational 
frame, so Ax'^ = 0 . Denoting the spatial part of the par¬ 
ticle 4 -velocity as n = 7/3, where /3 is 3 -velocity and 
7 = (1 — one obtains 


Ax* = A* Ax'^' (32) 

^ (1 + 7) (33) 

{Ax^f = AiAx'^AiAx'^ (34) 

~ ^jk T '^j'^k (33) 

(Ax)^ = (Ax')^ + (Ax' • u)‘^ (36) 

K{r) = No{2TTa‘^)~^/‘^exp 


To determine the normalization factor Vq, one can 
use that f (YYi^idxi) (detA)~^l^. The 

determinant (iet(A*A^) = det{Sjk + UjUk) = 7^, so 
J K(r)dAr = V0/7 = 1 and therefore Nq = 7. Finally 
we obtain 


K { r ) = {2 TTa ^) exp (33) 


One can check that this formula turns into the commonly 
used one if u = / 3 z)‘ Since 1 + — Iz 

in this case, we immediately retain the usual smearing 
kernel. Scalar products do not change under rotations, so 
the expression in the exponential does not change under 
rotations either. Rotating to the frame, where z axis 
is parallel to u one can see that our kernel is nothing 
else, but a proper Lorentz contraction in the direction of 
motion. 
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